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of  which  are  believed  to  be  previously  unsolved,  and  Green's  functions  are 
represented  using  the  special  functions  of  mathematical  physics.  A more 
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FOREWORD 


* In  this  thesis,  an  addition  formula  for  Green's 

I 

i 

i functions  of  linear  separable  elliptic  partial  dif- 

I 

I 

j ferential  equations  in  two  independent  variables  is 

derived.  The  formula  is  applied  to  generate  new  Green's 

! functions  from  known  Green's  functions. 

1 

This  thesis  is  in  manuscript  format;  it  consists 
of  one  manuscript  and  six  supporting  appendices.  The 
recommended  starting  point  for  reading  the  thesis  is  - 
Appendix  A:  Overview  of  the  Thesis. 
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AN  ADDITION  FORMULA  FOR 
GREEN'S  FUNCTIONS 


JEFFREY  S,  COHEN  and  JOHN  S.  PAPADAKIS 

Naval  Underwater  Systems  Center 

Nev  London  Laboratory 

New  London,  Connecticut  05320 


1 . INTRODUCTION 

One  of  the  most  general  approaches  to  solving  bound- 
ary value  problems  for  linear  differential  equations 
consists  of  expressing  the  solution  in  terms  of  an 

avixiliary  function  called  the  Green ' s function . The 

Green's  function  is  usually  defined  to  be  the  solution 

% 

of  the  adjoint  problem  with  a delta  function  as  a non- 
homogeneous  term.  The  particular  boundary  conditions 

affect  the  formulation  of  the  adjoint  problem  amd  thus, 
the  Green's  function  as  well.  The  Green's  function 
method  cam.  be  applied  to  elliptic,  parabolic,  amd  hyper- 
bolic linear  partial  differential  equations  as  well  as 
linear  ordinaury  differential  equations. 

This  paper  presents  a derivation  of  a formula 

that  gives  a form  of  the  Green's  function  of  a separ- 
aUsle  elliptic  partial  differential  equation  in  terms 

of  the  Green's  functions  for  two  .simpler  equations,  one 
of  which  is  elliptic  and  the  other  hyperbolic.  Because 


r 


of  the  presence  of  the  hyperbolic  equation  and  for  the 
saka  of  completeness,  we  shall,  in  section  2,  briefly 
review  the  definition  of  the  Green's  function  for  a 
Cauchy  problem  and  its  relationship  to  the  Riemann 
function.  Then,  the  Riemann  function,  which  is  known 
in  some  special  cases , can  be  used  to  find  the  Green ' s 
function  in  these  cases.  In  section  3,  the  formula 
will  be  formally  derived,  applied  to  the  reduced  wave 
equation,  and  be  further  justified.  In  section  4,  we 
shall  apply  the  formula  to  derive  what  we  believe  to 
be  a new  Green's  function. 

Before  proceeding  however,  it  is  interesting  to 
note  that  although  we  are  primarily  interested  in 
Green's  functions  for  elliptic  problems,  the  motivation 
' for  this  work  caune  from  recent  papers  on  hyperbolic 
problems.  A second-order  lineaur  hyperbolic  aquation 
in  two  independent  vauriables  may  be  solved  using  the 
Riemann  function.  In  1958,  Copson  [Ref.  1]  formally 
derived  a technique  giving  a possible  form  for  the 
Riejmamn  function  (which  Copson  called  the  Riemann-Green 
function) , assuming  the  variables  are  separable.  The 
technique  utilized  integral  transforms  which  generali- 
zed Riemann 's  own  use  of  Fourier  cosine  transforms.  In 
1964,  Mackie  [Ref.  2]  investigated  similar  problems  and 
presented  a close  relationship  between  the  Riememn 
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function  and  the  Green ' s function  for  the  Cauchy  prob- 
lem. Mackie  then  applied  Pourier-Beesel  trsinsforms  to 
find  the  Green's  function  for  Riemann's  original  example 
and  utilized  the  relationship  between  Riemann  and 
Green's  fiinctions  to  find  the  Riemann  functions.  Re- 
cently, Papadakis  and  Wood  [Ref.  3]  developed  an  addition 
formula  that  gives  the  Riemamn  function  of  a separable 
equation  in  terms  of  the  RiemcUin  functions  of  two  simpler 
equations.  In  deriving  the  formula,  they  employed  a 
refinement  of  Copson’s  integral  transform  method  to  find 
Green's  functions  and  then  utilized  Mackie 's  relation- 
ship to  find  the  Riemann  function.  In  this  paper,  the 
technique  is  further  developed  by  applying  it  to  Green's 
functions  for  elliptic  partial  differential  equations. 


2 . PRELIMINARIES 


a.  The  Green's  Function  for  Elliptic  Equations 
Consider  the  linear  elliptic  second-order  par- 
tial differential  equation  in  two  variables 

’ ’^xx  ^ ’^yy  2a(x,y)u^  - 2b(x,y)Uy 

+ c(x,y)u(x,y)  = f(x,y)  (2.1) 

with  suitable  boundary  conditions  on  a closed  curve  T 
(e.g. , u = 0 on  F) . The  Green's  function  for  (2.1), 
GCx»y*X,Y)  is  defined  as  the  solution  to  the  adjoint 
equation 

“ ^XX  + ®yy  “ 2(aG)^  + 2(bG)y  + cG 

- d-Cx-Xj  (TCy-Y)  (2.2) 

with  suitable  (adjoint)  conditions  on  the  boundary  F 
(e.g.,  G » 0 on  F) . If  D is  the  interior  of  F, 
the  solution  of  (2.1)  can  be  written  in  the  foirm 

u(X,y)  » I |G(x,y;X,Y) f (x,v)  dxdv  + [ B[u,G]  (2.3) 

j F 

where  B[u,G]  is  known  on  F (e.g.,  B[u,G]  *0  on  F) 

b.  The  Green's  function  for  Hyperbolic  Equations 
The  concept  of  the  Green's  function  is  also 

meaningful  for  the  Cauchy  problem 
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H[u]  * + 2a(X/Y)u^  - 2b(x,y)u^ 

+ c (x,y)u (x,y)  = f(x,y)  (2.4) 

with  initial  value  u and  given  on  the  line 

y » Vq.  As  in  the  elliptic  case,  the  Green's  function 
G(x,yyX,Y)  is  the  solution  of 

H*[G]  - ^xx  " ^yy  ” 2(bG)y  + cG 

» d{x-X)  (J(y-Y)  (2.5) 


subject  to  the  adjoint  boundary  conditions 

G - - 0 on  r (2.6) 

where  -2^  is  the  (outward)  normal  derivative  of  G on 
an 

r,  and  r is  a non-characteristic  curve  (i.e.,  the  ab- 
solute value  of  the  slope  of  r is  always  less  them  1) 

such  that  r and  a segment,  of  the  line  y * yg 
form  a closed  curve,  as  in  Pig.  1.  Prom  (2.5),  (2.6), 
cuid  the  eunount  of  freedom  we  have  in  choosing  T , it  can 
be  seen  that  G » 0 in  the  portion  of  the  interior  of 

outside  the  triangle  A,  with  vertices  at  points 
P,  A,  and  B.  Now  the  solution  u of  (2.4)  is  given  by 


u(X,Y) 


G(x,y;X,Y)  f (x,y) 


dxdy 


(u  G-uG  +2buG)  dx. 

y .y 


(2.7) 
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c.  The  Riemann  Function 

For  the  Cauchy  problem  described  above,  the 
Riemann  function  R(x,y;X,Y)  is  the  solution  of 

H*[R]  » 0 (2.8a) 

satisfying 


R^  + Ry.  » (a+b)R  on  y-Y  » x-X 

(2.8b) 

R - R “ (a-b)R  on  y-Y  = -(x-X) 

* y 

(2.8c) 

R(X,Y;X,Y)  » 1. 

(2.8d) 

Now  the  solution  of  (2.4)  is  given  by 

u(X,Y)  * - •^  [ ] R(x,y;X,Y)  f (x,y)dxdy 

+ j J (UyR-uRy+2buR)  dx 

» 

• 

+ i[u(A)R(A;P)  + u(B)R(B;P)]. 

(2.9) 

If  the  initial  values  of  u amd  Uy 

are  homo- 

geneous  (i.e.,  identically  0)  along  y = 

Yq , a compari- 

son  of  (2.7)  and  (2.9)  would  lead  to 

G(x,y;X,Y)  - - j R(x,y;X,Y) 

(2.10) 

inside  the  triangle  A.  If  the  initial  data  along 
y * 7q  is  nonhomogeneous , the  contour  integral  in  (2.7) 
requires  integration  of  terms  including  G and 
along  lines  across  which  G itself  is  discontinuous. 
This  generates  terms  analogous  to  those  inside  the 
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brackets  in  (2.9).  Thus,  (2.10)  holds  for  arbitrary 
initial  data.  (Except  for  a difference  in  notation 
between  a Green's  function  and  an  adjoint  Green's  func- 
tion, the  analysis  leading  to  (2.10)  is  due  to  Mackie 
[Ref.  2].) 
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3.  THE  ADDITION  FORMULA 
a.  Formal  Derivation 

Let  us  consider  the  separable  elliptic  equation 

* '^xx  ^yy  [c^(x)+C2(y)]u(x,y)  - f(x,y)  (3.1) 

and  attempt  to  find  the  "free-space"  Green's  function 
(i.e,,  assuming  c^^Cx)  and  Cj  (y)  are  analytic,  find 

a Green's  function  over  the  entire  (x,y)  plane).  Since 
L is  self-adjoint,  G(x,y;X,y)  must  satisfy 


L*[G]  = L[G]  » 5(x-X)  5(y-Y). 


(3.2) 


We  shall  attempt  to  find  an  expression  for  G in  terms 
of  the  Green's  function  G^  satisfying  the  hyperbolic 
equation 


Gi  - G^  + c (x)G  =*  5(x-X)  d(y-y) 

vv  -L  X 


(3.3) 


XX 


2uid  the  Green's  function  G2  satisfying  the  elliptic 
equation 


Gj  + G,  + Cj(y)Gj  - 5(x-X)  5(y-y) 
XX  ^ ^ 


(3.4) 

It  has  been  shown  by  PapadeUcis  and  Wood  [Ref.  3] 


that 


Rj^(x,y;X,Y)  = Rj^(x,y-Y;X,0)  = R^  (x,Y-y;X,0)  . (3.5) 

Recalling  (-2.10)  we  can  say  that  G^^  is  also  an  even 


, . .V.  - -K 
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function  of  y-Y  wherever  it  is  defined  for  both 
+ (y-Y)  amd  -(y-Y).  We  should  like  to  itiedce  a similar 
claim  for  G^.  If  we  assume  a condition  at  *,  call 
it  a "radiation  condition,"  ensuring  the  uniqueness  of 
G2*  then  we  can  claim  that 

G2(x,y;X,Y)  = Gj (x-X,y;0,Y)  H Gj (X-x,y;0,Y) . (3.6) 

Clearly,  any  form  of  (3.6)  satisfies  (3.4).  Thus,  if 

our  "radiation  condition"  is  also  satisfied  by  any  form 

of  (3.6),  the  claim  is  justified  by  the  uniqueness  of 

G2»  Since  we  are  proceeding  formally  at  this  stage,  let 

us  accept  these  hypotheses  and  (3.6)  as  est2U3lished, 

2 

with  this  final  observation.  If  ^2(7)  = k (the  case 
with  which  we  will  be  primarily  concerned  in  this  paper) 
auad  the  well-known  Sommer f eld  radiation  condition  [Ref. 
4]  is  assumed,  then  G2  is  unique  and  an  even  function 
of  x-X  as  desired. 

Now,  let  us  proceed  by  separating  variables  in  the 
homogeneous  form  of  (3.2),  leading  to  a pair  of  ordin- 
axy  differential  equations,  one  of  which  is 

9"(x)  + [Cj_(x)  + X^]9(x)  - 0 (3.7) 

2 

where  X is  the  sepeuration  constant.  Assume  there 
exists  a solution  9(x,X)  of  (3.7)  that  defines  a 
transform 
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f(X)  - J 9(x,X)F(x)  dx. 

Let  9(x,X)  be  the  inverse  transform  so  that 
F(x)  » I 9(x,X)f  (X)  dX. 

Applying  the  transform  to  (3.2),  integrating  by  parts 
twice  (ignoring  terms  evaluated  at  the  endpoints  of  the 
integral),  and  recalling  that  9(x,X)  is  a solution  of 
(3.7)  yields 

gyy  + [C2(y)  - X^]g  - 5(y-Y)  (3.9) 

where  g is  defined  by 

g(y;X,Y,X)9(X,X)  = | 9 (x,X) G(x,y;X,Y)  dx. 

Dtilizing  the  inverse  transform  gives 

G(x,y;X,Y)  » I 9(x,X)9(X,X)g(y;X,Y,X)  dX.  (3.10) 

Now  let  us  apply  the  same  procedure  to  (3.3). 
Separating  variables  leads  to  (3.7)  again  and  there- 
fore, the  same  solution  9(x,X).  Transforming  (3.3) 
yields 

gi  + X^g^_  - -6(y-Y)  (3.11) 

yy 

where 


g (y-Y;X,X)9(X,X)  = f 9 (x,X) G (x,y-Y;X, 0)  dx.  (3.12) 


r ' ' ■ 1 

11 

Recalling  boundary  condition  (2.6)  and  Rig.  1,  (3.11) 

can  be  solved  uniquely  with  t.he  initial  conditions 

g^^  '*■  0 and  g^^  -*-0  as  y ■*’ 

y 

g]_(y-Y;X,X)  - - gj-n  ^ (y-Y)  H(y-y) 

X 

I 

where  H is  the  Heaviside  step  function.  Inverse 
transforming  (3.12)  yields 

Gj^(x,y-Y;X,0)  -H(y-Y)  j d (x^)  d (X, \)  , (3.13)  | 

For  the  last  time,  separating  Veuriables  in  (3.4) 
leads  to  the  equation 

<J)"  (x)  + X^cp  (x)  » 0 

i 

with  a solution  of  the  form  e yielding  a Fourier 

transform.  As  before,  we  transform  (3.4)  and  find 

+ tc2(y)-X^]g2  - 5 (y-Y)  (3.14) 

where 

g2(y;X,Y,X)e  H j (x-X,y; 0 ,Y)  dx.  (3.15) 

—00 

! 

^ I Recalling  from  (3.6)  that  G2  is  an  even  function  of 

j x-X,  it  is  easily  deduced  that  g2  is  an  even  function 

I .fi  \ • f 

of  A and  independent  of  x.  Also  comp^u^ing  (3.14)  to 
(3.9)  leads  to  the  conclusion  that  gj  h g and  there- 
fore ' 


.V.  -A. 


L2 


t 


g(y;X,y,A)  = g(y,Y;X)  = g(y,Y;-\). 


Now,  inverse  transforming  (3.15)  and  using  the  fact 
that  G2  and  g are  even 

G- (x-X,y;0,Y)  » i cosX (x-X) g (v, Y; X ) dX.  (3.16) 

^ " J Q 

Let  us  formally  assume  that  the  integrals  in  (3.10) 
and  (3.13)  also  have  limits  from  0 to  *.  If,  in 
(3.13),  we  denote  y-Y  by  t and  investigate  for  t>0 


we  get 


G,  (x,t;X,0)  » -!  9 (x,X) 9 (X,X)cosXt  dX. 

Rewriting  (3.16)  with  t replacing  x-X 

If** 

* if’  j g(y»^;^)cosXt  dX. 

An  apolication  of  Parseval's  identity  to  (3.17) 
(3.18)  yields 


(3.17) 


(3.18) 

and 


G^  (x,t?X,0)G2(t,y;0,Y) 


dt 


9(x,X)  9(X,X)g(y,Y;X) 


dX 


1 

2 


G(x,y?X,Y) 


the  last  equality  coming  from  (3.10),  with  the  assumed 

limits  of  inteoration  from  0 to  «.  Thus  we  hope  to 
find  the  Green's  function  of  (3.1)  in  the  form 
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G(x,v;X,Y)  = -2  G (x, t;X, 0) G- (t ,y ; 0 , Y)  dt.  (3.19) 

^ 0 ■‘•t  ^ 

Equation  (3.19)  is  the  general  form  that  we  shall 
use  in  our  search  for  Green's  functions.  In  order  to 

get  a more  specific  formula,  one  must  investigate  the 
behavior  of  Gj_  and  and  the  effects  of  boundar- 

ies and  boundary  conditions  on  G.  For  the  time  being, 
we  shall  continue  on  with  our  assumptions  that  Cj_  (x) 

and  C2 (y)  are  well-behaved  functions  and  look  for  a 
"free-space"  formula.  In  section  4,  we  shall  analyze 
(3.19)  under  different  conditions. 

Now,  recalling  our  discussion  of  Green's  functions 
and  Riemann  functions  for  hyperbolic  equations,  and 
particularly  Fig.  1 and  Eq.  (2.10) 

I'O  for  y-Y  < |x-x! 

-|Rl(x,y;X,Y)  for  y-Y  ^ |x-x| 

where  the  Riemann  function,  R^,  is  defined  by  (2.3). 
Therefore,  we  can  write 

Gj_(x,t;X,0)  - -j  R^(x,t;X,0)H(t-|x-X|  ) 

where  H is  once  again  the  Heaviside  step  function 
Differentiating  G^_  with  respect  to  t; 


I 
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(x,t;X,0)  » -i.  R^(X,t;X,0)5  (t-|x-X|  ) 

- ^ R^^(x,t;X,0)H(t-|x-X|  ) . 


(3.20) 


Substituting  (3.20)  into  (3.19); 

G(x,y;X,Y)  » f {Rt5  -f  R,  H}G.,  dt 

J 0 ^ ^ 

m 

- R^(x,lx-Xl  ;X,0)G2(|x-x1  ,y;0,Y) 

00 

+ f Rt  (x,t;X,0)G-(t,y;0,Y)dt  . (3.21) 

Mx-X|  ^ 2 

Integrating  (2.8b&c)  with  a = b = 0,  as  in  this  case, 
we  get 

Rj^(x,  |x-xl  ;X,0)  - 1. 

From  (3.S) 


G2(|x-X|  ,y;0,Y)  =«  G2  (x-X,y;  0 ,Y)  . 

Lastly,  R^  and  G2/  as  they  appe2u:  in  the  integral  in 
(3.21),  are  even  functions  of  t so  that  the  integrand 

(with  its  t-derivative  term)  is  odd, and  we  can  drop 
the  absolute  value  from  the  lower  limit  of  integration. 
Recalling  (2.10)  we  can  now  write  (3.21)  as 


G(x,y;X,Y) 


G2 (x-X,y;0,Y) 


2 G,  (x,t?X,0)G2 (t,y;0,Y) 
^x-X  ^ 


(3.22) 

dt. 


Integrating  by  parts  auid  assuming  that  the  term  ^^^^2 
evaluated  at  t « « is  zero,  yields  the  alternate  form 
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G(x,y;X,y)  « 2 G,  (x, t ;X, 0 ) G,  (t,y;0,Y)  dt.  (3.23) 
^x-X  ^ 

Equations  (3.22)  auid  (3.23)  are  the  relationships  be- 
tween the  Green ! s functions  G,  Gj^,  and  G2  satisfying 
(3.2),  (3.3),  and  (3.4)  which  we  wished  to  derive  from 

our  general  form  (3.19).  Of  course,  all  our  work  so 
far  is  entirely  foinnal  and  thus,  justification  is  re- 
quired. Although  we  do  not  intend  to  produce  a rigor- 
ous proof  of  the  validity  of  the  formulas,  we  shall 
verify  that  they-  satisfy  some  of  the  properties  of 
Green's  functions  and  further  justify  some  of  our  as- 
sumptions. But  before  continuing  in  this  endeavor, 
let  us  convince  ourselves  that  it  is  worth  the  effort 
by  testing  the  formulas  on  a problem  for  which  we  know 
the  correct  answer. 

b.  A Test  Problem  - The  Reduced  Wave  Equation 
In  this  section,  we  shall  utilize  (3.23)  to 

2 

find  the  Green's  function  for  (3.1)  with  Cj_(x)  » -a 
2 

and  C2 (y)  = b , where  a and  b are  real,  positive 
constants  satisfying 

k^  * - a^  > 0.  (3.24) 

2 

Notice  that  we  have  broken  up  the  positive  constant  k 
as  the  difference  of  two  positive  consteuits  rather  than 

the  siun.  We  shall  comment  on  this  choice  later  in  the 
section. 
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Now  we  Ccm  rewrite  (3.2),  (3.3),  and  (3.4),  the 
equations  for  G,  G, , and  G^  - respectively, as 


+ I'  G * 5(x-X)  5(y-Y)  , (3.25) 

XX  yY 

G,  - G,  - a^G,  » 5 (x-X)  5 (y-Y) , (3.26) 

•‘•XX  •^yy 

and 

G + G + b^G  » 5(x-X)  5 (y-Y)  . (3.27) 

^xx  ^yy 

The  solutions  of  (3.26)  and  (3.27)  are  Jcnown  to  be 


(a/ (y-Y)  (x-X)  ^)  for  y-Y>lx-x| 


(x,y?X,Y) 


for  y-Y<|x-xl 


G^(x,yjX,Y)  » ^ (b/(x-X)^+(y-Y)^)  . (3.29) 

Substituting  (3.28)  and  (3.29)  into  (3.23)  and  then 
performing  the  required  differentiation  yields 

G(x,y;X,Y)  = Jq  (a/t^-  (x-X)  ^ (b/t^+  (y-Y)  ^ ) dt. 

x*X 

Jq  (bA^+  (y-Y)  ^)  • 

x-X 

bt  dt 


''t^t  (y-Y)  ^ 


Changing  the  integration  variable  from  t to 


3 » i^-(x-X)^ 
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yields 


00  y- — 

G(x,y;X,Y)  = ^ [ J.  (as)  (3.30) 

0 1 


where  r = / [x-X) (v-Y) ^ 


Utilizing  (3.24)  and  the  appropriate  integral  tables 

[Ref.  5,  pg.  358,  19.4  (3)  smd  Ref.  6,  pg.  706, 

6.596(6)1  yields 

G(x,y;X,Y)  » ^ (kr)  (3.31) 

the  well-lcnown  free-space  Green's  function  for  the  re- 
duced wave  equation  (3.25). 

Returning  to  the  remark  following  (3.24),  we  note 

2 2 *2 
that  if  k were  the  sum  of  a and  b . then,  in 

(3.28),  Gj^  would  consist  of  an  Ig,  a modified  Bessel 

function,  rather  than  a Jg.  If,  in  (3.30),  the  Jg 

term  were  replaced  by  ig,  the  integral  would  diverge 

since  Ig  increases  exponentially  for  large  arguments. 

Thus,  even  under  conditions  when  the  formula  is  appli- 

Cedsle,  each  problem  must  be  studied  with  great  care 

and  a certain  amount  of  ingenuity  may  be  required  in 

order  to  solve  the  problem. 

Let  us  note  that  the  formula  also  works  for  the 
2 2 

case  Cj^(x)  ■ -a  and  (y)  * -b  . The  analysis  is 

almost  identical  with  the  preceding  example,  with  the 
elliptic  Green's  function 


replaced  by 
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i-  K-,  the  modified  Hankel  function.  Reference  6 i 

IT  0 ! 

[pg.  706,  6.596(6)]  provides  the  required  integral. 

c.  Further  Justification  of  the  Addition  Formula 
In  order  to  verify  the  addition  formula  (3.19) 
or  formulas  derived  from  it,  such  as  (3.22)  and  (3.23),  j 

at  least  three  properties  of  our  representation  of  G ^ 

should  be  checked; 

> 

(1)  that  the  integrals  involved  converge; 

(2)  that  as  (x,y)  -*•  (X,Y)  , G possesses  the 
appropriate  (logarithmic'  singularity;  and 

(3)  that  for  (x,y)  (X,Y)  , G satisfies  the 

homogeneous  equation. 

Properties  (1)  and  (2)  should  be  verified  for  each  -< 

specific  problem  investigated;  that  is,  for  the  parti-  | 

cular  G,  and  G«  that  happen  to  arise.  For  for-  I 


y and  recalling  that 


G2  = -G2  - c^G^  (from  (3.4)) 

yy  XX 


G » G-  + 2 / G.  G2  dt  + 2c«f  G-  G^  dt.  (3.33) 
yy  2 ^x-X  ^t  ^tt  2Jx-x  It  2 


yy 

Adding  (3.32)  amd  (3.33) 


/ ^1  ^ 

^x-X  -^t 

r” 

+2  G,  G,  +G,  G,  )dt-2G~  (G,  ) 

^x-X  2tt  Itt  2t  2x  lx  y_Y»x-X 

00  00 

-c,G-C-(G.-2  G,  G,dt)  +2  (G,  G,  ) ^dt 

^ 2 2 2 It  2t  t 


“ 2G2  (G-t 

^X  •‘■X 


) . 


y-Y=x-X 


So,  using  (3.22) 


®XX  ®yy  (VC2)G 


2(Gi^G2^) 


- 2G-  (G, 

t»x-X  ^ 


) . 


y-Y»x-X 

Recalling  that  we  have  already  assumed  that  G^G2  ■*•  0 
as  t -►  « in  order  to  integrate  by  parts  ^ let  us  fur- 
ther assume  that  G,  G*  -►0  as  t •*■  <».  Then 

^t  2t 

G „ + ■*■  (C-+C_)G  ® — 2G0  (G.  +G.  ) » 0 

yy  12  2^  lx  ly  y.y.x-X 

by  C2.10)  and  integration  of  (2.8b,  with  a » b * 0) . 
Thus, (3.22)  and  (3.23)  satisfy  the  homogeneous  equation 
away  from  the  source  if 
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1 


lim  G, (x,t;X,0)G* (t,y;0,Y)  = 0 

t-^oo  -L  * 


and 


lim  Gi  (x,t;X,0)G2  (t,y;0,Y)  “ 
t-^oo  ^ ^ 


(3. 34) 


(3.35) 


These  conditions  were  met  by  the  Green's  functions  in 
our  test  problem. 

Let  us  note  that  in  some  cases  conditions  (3.34) 
and  (3.35)  may  be  unnecessarily  restrictive.  In  parti- 
cular (3.34)  was  used  to  integrate  (3.22)  by  parts  and, 
thus,  deriving  (3.23).  Then  we  used  the  equivalence  of 
(3.22)  amd  (3.23),  along  with  (3.35),  to  verify  that 
our  formula  did  indeed  yield  a solution  of  the  homo- 
geneous differential  equation.  However,  it  is  possible 
for  (3.22)  (or  (3.23),  for  that  matter)  to  be  valid 
while  (3.23)  (or  (3.22))  is  not. 

Consider  Laplace's  equation 


u, 


XX 


0. 


(3.36) 


With  c (x)  and  c (y)  both  0,  G.  and  G_  are  the 
1 2 12 

Green ' s functions  for 


u — u * 0 
XX  yy 


and 


u + u *0 
XX  yy 


respectively. 

Then, 


(3.37) 


(3.38) 
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-y  for  y-Y  ^ |x-x| 

G. (x,y;X,Y)  = ^ , (3.39) 

^ 0 for  y-Y  < Ix-Xl 

and,  letting  r = / (x-X) (y-Y) ^ 

* §7  (3.40) 

Substituting  in  (3.22)  we  find  the  well-lcnown  free- 
space  Green's  function  for  Laplace's  equation 

G(x,y;X,Y)  * ^ log  r ^ (3.41) 

since  G^^  ■ 0 inside  the  integral  appearing  in  (3.22). 
However,  (3.23)  diverges  for  the  same  situation. 

Now,  in  light  of  our  successful  general  verifica- 
tion of  property  (3)  and  our  successful  applications  of 
the  formulas  in  the  test  example  and  Laplace ' s equa- 
tion, let  us  take  a heuristic  look  at  properties  (1) 
amd  (2)  by  placing  certain  restrictions  on  c^^(x)  and 
C2 (y)  and  the  boundary  conditions  associated  with  (3.1) 
For  the  formal  derivation, we  assumed  that  c^  and  C2 
were  analytic  in  the  entire  plane.  Now, let  us  modify 
this  assumption  so  that  rather  than  the  whole  plame, 
we  aure  looking  at  an  infinite  domain  with  boundaries 
separable  in  x amd  y (e.g.,  the  half-plane  x > 0). 
Also  let  us  assume  that  any  singularities  in  c^  amd 
C2  appear  outside  the  domain  of  interest.  Finally, 
assume  that  as  x and  y approach  ®,  c^^(x)  and 
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C2(y)  approach  real  constant  values.  Without  loss  of 
generality,  it  can  be  assumed  that  c^^  approaches  a 
negative  constant  since  we  can  subtract  a constamt  from 
c^  and  add  it  to  C2  and  not  alter  the  total  coeffi- 
cient in  (3. ] ) . 

With  these  assumptions  on  c^  and  , we  can 
hope  that  and  G2,  satisfying  (3.3)  and  (3.4), 

will  behave,  near  *,  like  and  (or  Kq), 

the  Green's  functions  satisfying  (3.26)  and  (3.27). 

From  the  test  problem,  we  know  that  the  integral  (3.30) 
converges  and  thus  hope  that  the  integrals  .containing 
Gj^  amd  G2/  (3.22)  amd  (3.23),  will  also  converge  and 
that  the  integration  by  parts  leading  from  one  form  to 
the  other  will  be  justified.  If  necessary,  further 
restrictions  on  how  fast  c^  and  C2  approach  their 
limiting  vadues  could  be  made  to  insure  asymptotic 
knowledge  of  G^^  amd  G2,  or  perhaps  entirely  differ- 
ent types  of  assumptions  are  necessary.  However,  since 
the  derivation  of  the  addition  formula  has  been  formal, 
we  shall  leave  the  problem  of  convergence  at  this  point. 

We  shall  cope  with  the  problem  of  the  singularity 
as  (x,y)  -►  (X,Y)  in  much  the  same  vein.  Since  we 
have  assumed  that  c^  and  C2  are  well-behaved,  we 
shall  also  assume  that  they  can  be  treated  locally 
as  constants.  With  these  assumptions  we  can  hope  that 
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counterparts  in  the  test  problem  cuid  exhibit  the  ap- 
propriate singularity. 

While  these  arguments  indicate  that  properties 
(1)  and  (2)  will  be  satisfied  by  certain  types  of  in- 
teresting problems,  let  us  recall  that  the&e  proper- 
ties should  be  checked  for  each  individual  problem. 
Property  (3) , however,  has  been  directly  verified 
under  general  conditions. 
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4.  THE  GREEN'S  FUNCTION  FOR  + u + [k^-  ^.q 

XX  yy  ^2  * 

a.  Derivation  of  the  Green's  Function  using  the 
Addition  Formula 

In  this  section, we  shall  finally  apply  the  addi- 
tion formula  to  a non-trivial  example.  To  the  best  of 

our  knowledge,  the  Green's  function  found  in  this  sec- 
tion has  not  been  previously  presented. 

Let  us  consider  the  equation 


L[u] 


u, 


XX 


+ u [k^  - 

x2 


yy 


(4.1)  I 

1 1 


on  the  half-plane  x > 0,  where  m is  a positive  inte- 
ger. We  shall  attempt  to  find  the  Green's  function 
G(x,y;X,y)  for  (4.1)  on  the  half-plane  x > (J  satis- 
fying 


L[q]  - 

5(x-X)  5(y-y) 

(4.2) 

G » 0 

on  X » 0 

(4.2a) 

G 0 

as  1 (x,y)  1 « 

(4.2b) 

As  (4.1)  seems  to  suggest  we  shall  apply  our  for- 
mula with  c,  (x)  * - ” and  c- (v)  » k^.  Thus,  we 
must  find  the  Green's  function  G^  for 


u - u 
XX  yy 


m(m-fl) 


As  we  have  seen  previously,  G2 


(4.3) 

must  satisfy  (3.25);  so 


i . 
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G2(x-X,y;0,Y)  - ^ (kr) 
where  r = / (x-X) (v-Y) ^ . 


In  order  to  avoid  an  unnecessary  digression,  we 
shall  show  in  Appendix  B that 


^ -i  P (?)  for  I x-x| <y-Y<x+X 


(x,y-Y;X,0) 


2 m 
0 


elsewhere 


(4.4) 


2 2 2 

where  ? = ^ — and  P_  is  the  Legendre  oolv- 

2xX  ™ ^ . 

nomial  of  the  first  kind  (see  Pig.  2) . 

Now  that  and  G2  are  known,  we  can  return  to 

our  general  form  of  the  addition  formula,  Eq.  (3.19). 
Then  (4.4)  may  also  be  written  in  the  form 

G,(x,t;X,0)  = -1  ) H ( t~  I x-X  I ) H(x+X-t). 

1 2m  2xX 

Noting  that  P (?)  is  the  Riemann  function  R. (x,y;X,Y) 
m,  1 

for  (4.3),  we  can  also  write 


Gj_(x,t;X,0)  » R^(x,t;X,0)  H(t-lx-X|)  H(x+X-t).  (4.5) 


Differentiating  with  respect  to  t 

Gi^(x,t;X,0)  » -j  Rj_^  (x,t;X,C)H(t-|x-Xl  )H(x+X-t) 

R^(x,t;X,0)6 (t-|x-x| )H(x+X-t) 
+j  R^(x,t;X,0)H(t-|x-X| )5 (x+X-t) . 


Substitution  into  (3.19)  yields 
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G2<it 


G2<it 


I H (x+X-  I x-X  I ) G2 

1 y-Y» I x-X I y-Y« | x-X | 


V-Y-X+X 


H(x+X-|x-X| )G2 


I y-Y»x+X 


Since  X emd  X are  both  positive,  the  Heaviside 


functions  are  both  1.  Since  is  the  Riemann  func- 


tion for  (4.3),  R^  « 1 along  the  characteristic  line 


y-Y  » x-X;  a fact  we  could  also  have  found  by  noting 


that  P., 


y-Y«  x-X 


? (1)  « 1.  Although  we  also  know 


y-Y-x+X 


P_(-l)  * (-1)“*,  we  shall  not  use  this  fact 


at  this  stage  of  the  analysis  because  it  appears  to  be 


too  specialized  (i.e.,  cannot  be  deduced  from  proper- 


ties of  the  Riemann  function]  to  include  in  a general 


form  of  the  formula.  How,  since  the  integrations  over 


the  discontinuities  of  G^  have  been  t2dcen  into  ac- 


count, we  can  replace  R^^  by  “2Gj^  and  using  the  now 

faunili2u:  argument  that  G2  is  an  even  function  of  t 
and  that  the  integrand  in  (4.6)  is  odd,  we  get 


.x+X 

G(x,y;X,Y)  » -2j.  G-i  (x,t;X,0)G- (t,y;0,Y)dt 
J x-X  2 


+ G2(x-X,y;0,Y) 


+ 2Gj^(x,x+X;X,0)G2  (x+X,y;0,Y)  . (4.7) 


Equation  (4.7)  is  the  counterpart  of  (3.22).  As 
a matter  of  fact,  in  this  case  we  did  not  need  to  go 

all  the  way  back  to  (3.19).  We  could  have  substituted 
our  expression  for  ^l^'  (4.5),  into  . (3. 22)  and  gotten 

(4.7).  However,  because  we  have  not  proven  a theorem 
as  to  the  apoliczdaility  of  the  formula,  we  gain  more 

confidence  by  reanalyzing  (3.19)  them  by  performing 
blind  substitutions.  Now,  integrating  (4.7)  by  parts 

yields  the  counterpart  of  (3.23) 

.x+X 

G(x,v;X,Y)  -2  G.  (x,t;X,0)G2^(t:,y;0,Y)dt.  (4.8) 

J x-X  ^ 

Direct  differentiation  of  (4.7)  and  (4.8),  in  a 

manner  analogous  to  that  performed  in  section  3c,  veri- 
fies that  G is  indeed  a solution  of  (4.1)  for 
(x,y)  ^ (X,Y).  Let  us  note  that  because  we  have  finite 
limits  of  integration,  we  no  longer  need  conditions 

(3.34)  and  (3.35). 


Substituting  our  known  forms  of  G^  amd  G2  into 
(4.7)  and  (4.8),  with  the  added  notational  convenience 
r = / (x+X) (y-Y) ^ , yields  ’ 
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I 


G(x,v;X,Y)  * f P (kX2+(y-Y)^)dt. 

4i  m 2xX  Ot 

(4.10) 

So,  for  any  positive  integer  m,  we  have  shown' 
that  the  addition  formula  yields  the  convergent  inte- 
grals (4.9)  and  (4.10)  which  also  satisfy  (4.1)  for 

•1) 

(x,y)  ^ (X,Y).  If  we  now  differentiate  the  Hq 

term  as  indicated  in  (4.10)  and  then  perform  the  change 

of  variables 

• 3 - (y-Y)2 

we  get 

G(x,v;X,Y)  « ^ p (ks)kds. 

4i  m 2xX  ^ 

(4.11) 

In  this  form,  we  shall  show  that  G has  the  appropriate 
singularity,  G 0 as  x 0,  and  G -►  0 as  r -►  *>. 

b.  Verification  of  the  Singularity  of  the  Green's 
Function 

For  notational  convenience  in  investigating 
the  behavior  of  G neeu:  the  singular  point  r * 0, 
let  us  write 

j _ X ■►X  -Ky-y)  -s  , r -cr  -2s  _ 

2xX  r^-r^ 

the  identity  following  directly  from  the  definitions 
of  r and  r.  From  (4.12)  it  is  evident  that,  for 


j 

i 


j 


j 

i 


1 

\ 
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0 < r < s and  for  s -*•  0/  ^ -^1.  Therefore 


lim  P_(C) 
0<r<3^0 


in 


1 


by  the  continuity  of  all  the  functions  involved.  So 
given  e > 0 arbitrarily  small,  choose  5 > 0 such 


that 


0 < r < s < 5 — ^ 1-e  < P^(C)  < 1 (4.13) 

and  such  that 
Yj^(ks)  < 0 

is  simultaneously  satisfied,  where  of  course  is 

the  first  order  Bessel  function  of  the  second  kind  (i.e.. 

Now  for  0 < r < 5 rewrite  (4.11)  as 
G » I Pj^(c)J^(ks)kds  + J 1*^  Pjjj(!;)Y^(ks)kds 

+ j I P (?)Y, (ks)kds.  (4.14) 

4 j 5 m i 

The  first  and  third  integrals  in  (4.14)  cure  well-behaved 
and  approach  constant  values  as  r -»•  0.  Recalling 

(4.13),  we  can  apply  the  mean  value  theorem  to  the  re- 
maining integral  resulting  in 

G * [^^(kr)  - YQ(k5)]  + constant 
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where  C*  is  evaluated  at  s*  e (r,6).  We  know  that 
for  any  s e (0,5)  , 1-e  ^ Pjjj  ^ 1 so  that 

j(l-£)YQ(kr)  + constant  ^ G £ j Yg (kr)  + constant. 

Thus,  as  r -*■  0,  G behaves  like  j Yg  (kr)  or  equival- 
ently ^ log  r (Ref.  7,  pg.  360,  9.8.1],  as  desired. 

We  have  now  verified  conditions  (1) , (2) , and  (3)  of 
section  (3c),  that  is,  the  integrals  in  our  representa- 
tion converge,  that  G possesses  the  correct  (loga- 
rithmic) singularity,  and  that  except  at  the  singularity 
G satisfies  the  homogeneous  partial  differential  equa- 
tion. All  that  remains  is  an  investigation  of  "the  boun- 
daury  conditions. 

c.  Verification  of  the  Boundary  Condition  along 

X * 0 

Approaching  the  singular  line  x = 0 in  the 
half-plane  x > 0,  we  have  0 < r < r.  In  this  case, 
for  r ^ 3 <_  r,  as  the  limits  of  integration  in  (4.11) 

indicate,  we  observe  from  (4.12)  that  "1  <,  ? <.  1 «ind 
therefore  <_  1 [Ref.  8].  We  can  now  write 

as 

< lim  7 f (ks)  Ikds  - 0 , ^ 

“ x-O  4 J r ' 1 ' 

since  as  x-^0,  r-*-r  and  r is  strictly  positive 
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(as  a matter  of  fact  r^X>0  on  x=0)  and, thus, 
|H^^^(ks)|  is  a continuous  bounded  function  in  the 
interval  from  r to  r. 


d.  Verification  of  the  Condition  at  Infinity 
, Lastly,  we  investigate  the  behavior  of  G as 
r ->■  ®.  Once  again,  0 < r < r and  I — ^ 

that,  from  (4.11), 


<3|  <j\  (ks)  |kds 

*'  r 

_<  (kr)  I (r-r)  . 


(4.15) 


since  |H^  (kr) 1 is  a monotonically  decreasing  func- 
tion. (Ref.  6,  pg.  969,  8.478],  Prom  the  definitions 
of  r and  f we  find  that 

0 < J-r  » r[(l  + l2Sjl/2  _ ^ 


For  r sufficiently  large,  simple  meuiipulations  show 
that 


0 < r-r  < 2X  + O(-) 
— r 


(4.16) 


For  r l2u:ge,  it  is  also  known  [Ref.  7,  pg.  365, 
9.2.28]  that 


(kr) I - 0(r“^/^)  . 


(4.17) 


Substituting  (4.16)  euid  (4.17)  into  (4.15)  yields 


....  «,  * "i  . ■ %.  * •» . .'*?•  - - • ' 4 


|g|  < 


2md,  thus , G -*•  0 as  r -*•  ®.  Hence,  we  have  demonstra- 
ted that  (4.11)  satisfies  (4.2),  (4.2a),  and  (4.2b). 

1 

e.  A Closed  Form  for  m = 1 

For  the  case  m » 1,  (4.11)  becomes  easy  to 
integrate,  utilizing  the  facts  that  P^(C)  » 5 and 

J (ks)kds  = (ks)  » |2.  (ks) 

- (ks)  (4.18) 

[see  Ref.  6,  pg.  634,  5.52(1)  and  pg.  967,  8.471]. 

From  the  first  equality  in  (4.18)  we  integrate  (4.11) 
to 


G(x,y;X,Y) 


^ 2£_±2L±IZzIi_[H^^^  (kr)-Hjj^^  (kJ)  ] 
2xX 


+ IT  (kr)-r^H^^^  (kJ)  ] 


(4.19) 


+ — (kr)-r^H^^^  (kf)  ] . (4.20 

4i  r^-r2  ^ 2 


Utilizing  the  second  form  of  (4.18)  yields 


G(x,y;X,Y)  - ^ (kr)  + ^ H^^^  (kr) 


^ [rH^^^kr)-rH|^^  (kr)  ] . (4.21) 


k (r2-r2) 
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f.  The  Case  k = 0 

Using  standard  ascending  series  forms  for 

J (z)  2uid  Y (z)  [Ref.  1,  og.  360,  9.1.10  and  9.1.11]  , 
n n 

it  is  easy  to  show  that 

(z)  » — log  z + — + Y~loq2]  + 0(z^log  z) 

U TT  IT 


z 

Substituting  these  forms  into  (4.19)  amd  letting  k -»■  0 
yields 

G(x,y;X,Y)  » 1 (y-jQ..) log  r/J  + i-  . (4.22) 

2TT  2xX  ^ 

Thus,  (4. 22)  is  the  Green's  function  for 


u + u - — u » 0 
MC  yy  ^2 


(4.23) 


satisfying  the  boundary  conditions  G -*•  0 as  x -*>  0 


and  as  [ (x-X,y-Y) 


».  It  is  interesting  to  note 


that  (4.22)  displays  precisely  the  same  form  as  a 

fundamental  solution  as  defined  by  Garabedian  [Ref.  9] . 
If  we  allow  k to  approach  0 through  positive 


values  in  (4.11),  we  get 


(4.24) 


(4,24)  can  also  be  derived  directly  from  the  addition 
formula  (4.8)  where  Gj^  still  is  the  Green's  function 
for  (4.3),  but  Go  « ^ log  r is  the  Green's  function 

A 2^ 


for  Laplace's  equation  (3.36).  In  the  case  m » 1, 
(4.24)  integrates  directly  to  our  solution  (4.22). 
Since  is  simply  a polynomial,  a straightforward 

integration  of  (4.24)  may  be  performed  for  all  posi- 

tive  integers  m,  yielding  closed  form  Green's  fvmc- 
tions  for 


+ u 


m(m+l)  .. 


yy  X* 

which  satisfy  the  boundary  conditions  G 


as 


X -►  0 as  r ®. 


V 


Show. 
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5.  SUMMART 

In  section  4,  we  applied  the  addition  formula  to 
the  equation 

u^^  + u + [k^  - u = 0 (5.1) 

Xvc  X" 

and  verified  that  the  integral  expression  we  derived 
was  indeed  a Green's  function.  For  the  particular  case 
ms  1,  we  were  aJsle  to  express  the  solution  in  a closed 
form.  In  addition,  we  succeeded  in  letting  k -*•  0 in 
the  general  integral  form  yielding  Green's  functions 
for  (5.1)  with  k * CT.  While  no  claims  have  been  made 
about  the  uniqueness  of  our  solutions,  we  feel  that, 
for  these  equations  the  presentation  of  solutions,  in 
such  relatively  simple  forms,  possessing  the  correct 
singularity  and  homogeneous  boundary  conditions  is, 

in  and  of  itself,  of  great  interest  and  importance. 

Furthermore,  the  addition  formula  derived  in  sec- 
tion 3 yields  a new  method  of  searching  for  Green's 
functions  for  separable  elliptic  partial  differential 

equations.  Of  course,  some  analysis  must  be  performed 
for  each  problem  encountered,  but  ovir  successful  ap- 
plication of  the  formula  to  the  'equations  studied  in 
sections  3 and  4 encourages  us  to  attempt  to  apply 
the  formula  to  other  cases. 
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Appendix  A:  OVERVIEW  OF  THE  THESIS 

1 . INTRODUCTION 

The  major  point  of  the  thesis  is  the  presentation 
of  a technique  which  yields  formulae  for  Green's  func- 
tions of  linear  separable  elliptic  partial  differen- 
tial equations  in  two  independent  variables.  The  for- 
mulae require  knowledge  of  the  Green's  function  (or 
Riemann  function)  of  a hyperbolic  equation  and  the 
Green's  function  of  an  elliptic  equation.  Both  the 
hyperbolic  amd  elliptic  equations  utilized  in  the  pro- 
cess are  substantially  simpler  than  the  original  equa- 
tion. Thus,  known  Riemann  and  Green's  functions  may 
be  used  to  generate  formulae  for  new  Green's  functions. 
Once  a formula  for  the  Green's  fionction  is  known,  one 
can  usually  justify  the  result  rigorously  using  stan- 
dard analytical  techniques.  In  addition,  the  overall 
process  is  applied  in  the  thesis  to  find  specific 
Green's  fui*ctions  satisfying  homogeneous  boundary  con- 
ditions which  we  believe  to  be  new  solutions . 
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2.  CONTEXT  OF  THE  THESIS 

The  memuscript  deals  with  the  self-adjoint  equa- 
tion 

'^xx  C2(y)]u  = 0.  (A.l) 

In  order  to  find  the  Green's  fxinction  G(x,y;X,Y) 
for  (A.l)^we  investigate  the  Green's  function 
Gj^(x,t;X,0)  of  the  hyperbolic  equation 

u^  - U(.|.  + Cj^(x)u  » 0 (A.  2) 

amd  the  Green's  function  G2(t,y;0/Y)  of  the  elliptic 
equation 

“tt  “yy 

The  motivation  for  calling  our  result  an  addition  for- 
mula comes  from  the  fact  that  formally  adding  (A. 2) 
and  (A. 3)  yields  (A.l). 

Employing  separation  of  variables,  tramsform 
techniques  for  solving  ordinary  differential  equa- 
tions, and  Parseval's  identity  from  Fourier  transform 
theory,  we  derive  the  addition  formula 

G(x,y;X,Y)  - G2(x-X,y;0,Y) 

OO 

-2  f G,  (x,t;X,0)G_ (t,y;0,Y)dt. 

' x-X  t ^ 


(A.  4) 
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Integrating  by  parts  cind  assuming  that  the  term  Gj_G2 

evaluated  at  t » ® is  zero,  yields  the  alternative 
form: 

G(x,y;X,y)  » 2[  G,  (x,t;X, 0) G,  (t,y; 0, Y) dt.  (A. 5) 

The  significance  of  formulas  (A. 4)  and  (A. 5)  is 
apparent.  If  c^^  (x)  is  such  that  the  Riemann  func- 
tion of  (A.  2)  is  known  or  can  be  found  and  ^2^.7)  is 
such  that  the  Green's  function  of  (A. 3)  is  known  or 
can  be  found,  one  cam  substitute  them  into  (A. 4)  and 
(A.  5)  and  then  directly  verify  that  G satisfies  the 
requirements  of  a Green's  function.  As  an  example,  we 

let  c-t  (x)  « - ^ so  that  G,  is  known  to  be 

^2  * 

an  m-th  degree  Legendre  polynomial  of  the  first  kind 

2 

(see  Appendix  B for  the  precise  form)  and  C2 (y)  » k , 

so  that  G2  is  a 0-th  order  Hankel  function  of  the 
first  kind,  (kr) . Substituting  these  into  (A. 5) 

we  get  am  expression  (4.11  in  the  manuscript)  for 
the  Green's  function  G(x,y;X,Y)  of  the  equation 

’^xx  '^yy  ] u = 0.  (A. 6) 

In  this  form,  we  prove  that  G is  a solution  having 
the  correct  (logarithmic)  singularity  and  satisfying 

homogeneous  conditions  along  the  singular  line  x « 0 
and  varnishing  at  infinity.  We  believe  this  result  to 


? 
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be  new,  at  least  in  the  particularly  simple  forms  pre- 
sented in  the  manuscript.  In  addition,  we  perform  the 
necessary  integrations  and  get  closed  forms  for  the 

Green's  function  of  (A. 6)  for  the  case  m » 1 and  k 
arbitrary  amd  for  the  case  k » 0 and  m arbitrary. 

At  this  point  in  the  thesis  development,  it  was 
felt  that  the  results  should  be  presented  in  the  open 
literature  for  the  use  of  other  mathematiciams  work- 
ing in  this  field.  However,  two  relatively  straight- 

forwaurd  questions  are  ignored  in  the  manuscript.  These 
deficiencies  are  eliminated  in  Appendices  C and  D.  In 
Appendix  C the  results  of  the  mamuscript  aure  generalized 
from  self-adjoint  equations  in  the  form  (A.l)  to  non- 
self-adjoint equations  in  the  form 

'^xx  '^yy  ■*'  2a(x)u^  + 2b(y)Uy  + [c^^  (x) +C2  (y)  1 u»0  (A.  7) 

amd  some  non-self-adjoint  problems  are  solved.  In 
Appendix  D we  solve  a problem  in  which  (y)  is  no 
longer  constant.  As  a matter  of  fact  we  use  the  results 
of  the  manuscript  to  find  the  Green's  function  for 


XX 


+ U - + Sl2±iLi  u 


yy 


(A.8) 


X**  y“ 

in  the  first  quadrauit  of  the  plame  satisfying  homogene- 
ous conditions  along  the  coordinate  auces.  Once  again, 
this  is  believed  to  be  a new  result.  The  fact  that  we 
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could  solve  otherwise  intractable  problems  like  these 
(i.e.,  Eqs.  (A. 6)  eind  (A. 8))  indicate  the  power  and 
utility  of  the  addition  formula. 
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3.  SPECULATIVE  DISCUSSION 

Having  derived  an  addition  formula  for  Green's 
functions  of  elliptic  equations,  one  begins  to  wonder 

how  the  result  can  be  generalized.  In  terms  of  func- 
tional analysis,  the  partial  differential  equation 

is  a linear  operator  and  the  process  of  finding  the 
Green's  function  basically  consists  of  inverting  the 
operator.  How  cam  we  characterize  those  linear  opera- 
tors which  admit  an  addition  formula  solution  (to 

itself  or  to  the  inverse  problem) ? At  the  present 
stage  of  development  we  cannot  answer  this  question. 

In  order  to  get  a hint  towards  a direction  to  pro- 
ceed, we  would  like  to  find  other  operators  for  which 
addition  formulae  may  be  found.  An  addition  formula 
for  Riemamn  functions  of  linear  separable  hyperbolic 

equations  in  two  independent  variables  has  been  found 
by  Papadakis  amd  Wood  [Reference  '3]  and  their  results 
are  stated  in  Appendix  E.  However,  we  still  do  not 
have  enough  examples  to  guide  us. 

In  section  3c  of  the  mamuscript,  it  is  proven 
that  the  addition  formulas  (A. 4)  and  (A. 5)  are  indeed 

solutions  of  (A.l)  away  from  the  singularity.  Follow- 
ing this  amadysis,  it  would  appeau:  that  if,  instead  of 

using  the  Green's  function  G2  of  (A. 3),  we  used  any 
solution  to  the  homogeneous  equation,  we  get  a solu- 
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tion  to  the  homogeneous  equation  (A. 4).  Thus,  we  might 

I 

be  able  to  develop  an  addition  formula  for  solutions  i 

rather  them  for  Green's  functions.  Now,  totally  enter- 
ing the  realm  of  speculation,  it  is  felt  that  cm  addi- 
tion formula  could  be  derived  for  elliptic  equations  in 
3 dimensions  in  terms  of  a hyperbolic  and  elliptic  } 

decomposition,  or  possibly  in  terms  of  a decomposition 
consisting  of  two  parabolic  equations.  It  is  also  ^ 

guessed  that  certain  types  of  non-separadsle  equations 
might  admit  some  form  of  an  addition  formula.  Further 
research  into  these  specific  areas  cmd  into  functional 

I; 

analysis  techniques  (such  as  those  utilized  in  referen- 
ces 4 and  5)  is  certainly  an  interesting  avenue  of 
study.. 
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Appendix  3:  THE  GREEN'S  FUNCTION  FOR 

m (m+1) 


u - u 
XX  yy 


u 


Consider 


L[G]  - - 2^2^  G » 5(x-X)  5 (y-Y)  . (B.l) 

XX  yy  ^2 

1/2 

Multiply  both  sides  of  (B.l)  by  x ' *^nH*l/2^^*^ 
integrate  with  respect  to  x from  0 to  •.  Aftor 
some  manipulation  (including  two  inteorations  by  parts 

and  an  application  of  Bessel's  ordinaury  diffsrsntisl 
equation ) 


where 


.00 

g(y;X,Y,X)  » J G (x,y;X,Y)  dx.  (B.3) 


Recalling  our  previous  discussions  on  the  behavior  of 

Green's  functions  for  hyperbolic  equations  and  the  simi- 
lar form  of  (3.11),  the  solution  of  (B.2)  is 

g(y-Y;X,X)  - <^^)  H (y-Y)  . (B.4) 

Noting  from  (B.3)  that  g is  the  Fourier-Bessel  trcuis- 
-1/2 

form  of  X G,  we  can  inverse  tr2msform  (B.4)  and 
1/2 

multiply  by  X ' to  yield 


G(x,y;X,Y) 


- -(xX)^^^H 


(y-Y)f  sinX  (y-Y)J_.l(Xx)J  i(XX)  dX. 
J Q m+^  m+y 
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Reference  to  an  appropriate  table  of  integrals  [Ref.  6, 
pg.  732,  6.672  (1)]  yields  the  desired  result,  which 
is  Eq.  (4.4)  for  m a positive  integer.  For  m not 
an  integer. 
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Appendix  C:  NON- SELF-ADJOINT  EQUATIONS 

Although  the  memuscript  deals  only  with  linear 
elliptic  partial  differential  equations  which  are 

separable  and  self-adjoint,  the  non-self- ad joint  case 
may  be  reduced  to  self-adjoint  form  by  a transforma- 
tion of  the  dependent  vzuriable.  Let  us  consider  the 
more  general  problem 

* Gj^+Gyy+2a  (x)G^+2b  (y)  G^+[Cj^  (x) +Cj  (y)  ]G(x,y;X,Y) 

* 5 (x-X)  6 (y-Y) . (C.l) 


Let  G(x,y;X,Y)  be  defined  by 

G(x,y;X,Y)  ■ exp{-f  a(t)dt-f“b(t)dt}G(x,y;X,Y)  . (C,2) 

^ Y 


Under  this  transformation  of  variables  , L 
trems formed  operator  L by 

L[Gl  = exp{-[  a(t)dt  - b(t)dt>L[G] 
^X  Y 


where 


defines  the 


(C,3) 


L[G]  » G^^  + G + tc,  (x)  + c-(y)]G.  (C.4) 

XX  yy  X z 

Eq.  (C.4)  is  precisely  in  the  form  dealt  with  in 
' the  manuscript.  Now  if  G satisfies 

L[G]  - 5(x-X)  5(y-Y)  (C.5) 

then  Eq.  (C.l)  is  satisfied  since  the  exponential  term 


in  Eq.  (C.3)  becomes  unity  when  (x,y)  » (X,Y). 
Returning  to  Eq.  (C.4)  we  note  that 


c (x)  » c,  (x) 
1 1 

and 

£2(7)  “ 


[a(x) - a'  (x) , 
[b(y)l^  - b*  (y)  . 


(C.6a) 


(C.6b) 


Let  L have  the  form  of  the  exeunple  in  section  4 

of  the  memuscript,  that  is 

LIG]  - + Gyy  + [k^  - ]G  , (C.7) 

Eq.  (4.11)  of  the  manuscript  gives  us  the  expression 

for  the  Green' s function  G satisfying  G » 0 on 
^ 2 2 

X » 0 and  G -►  0 as  x +y  “. 

Hence  we  can  generate  Green ' s f iinctions  for  the 

following  non-self -ad joint  equations; 

L.  [G.  1 * G^  + Gj^  — — ~ G]_  + k (C.8) 

■*"  XX  yy  * X 

amd 

2m 

L2[G2]  * G2  t G2  “ G2  hG2  • (C.9) 

^ XX  yy  X y 

The  Green's  function  for  (C.8)  is 


Gj^(x,y;X,Y)  « (|)®  G(x,y;X,Y) 


(C.IO) 


and  for  (C.9}  is 

G2(x,y;X,Y)  - (|)®  e“^^^“^^G(x,y;X,Y)  . (C.ll) 

Gj^  and  62  are  both  0 on  the  singular  line  x « 0 
but  they  no  longer  varnish  at  infinity. 
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Appendix  D:  THE  GREEN'S  FUNCTION  FOR 


u +u  - . 0. 


XX  yy 


‘In  section  3 of  the  manuscript/  the  addition  for- 
mula vras  derived  for  equations  of  the  form 


'^xx  'Vy  [c^(x)+C2(y)]u  = 0. 

However,  in  all  the  examples  in  the  manuscript,  C2 
was  identically  constant.  In  this  appendix  we  shall 
find  the  Green's  function  G(x,y;X,Y)  defined  in  the 
first  quadrant  of  the  x,y  pleme  (x  > 0,  y > 0) 
satisfying 

^xx  + “ [5Li2^  + n (ntllj  g , 6 (x-X)  5 (y-Y)  (D.l) 

yy  x^  - y2 

and  subject  to  the  boundary  conditions 


G ■ 0 on  X » 0 amd  on  y * 0. 


Following  the  analysis  of  section  4 of  the  mamu- 
script,  we  consider  the  Green's  function  Gj^(x,t;X,0) 

of  the  hyperbolic  equation 


“xx  - “tt  - 2^  “ - ° 
which  we  found  to  be  (see  Appendix  B ) 


(D.2) 


G^  (x,t;X/0) 


.ip 

2 m 2xX 


for  x-X  <t  < xtX 


(D.3) 


52 


and,  for  m a positive  integer,  0 elsewhere. 

Next  we  need  the  Green's  function  G2(t,y;0,Y) 
for  the  elliptic  equation 


U,,  * - 0. 


(D.4) 


The  solution  to  this  problem  with  boundary  condition 


G-  » 0 on  y « 0 


(D.5) 


was  the  subject  of  section  4f  of  the  manuscript  and 


with  a change  of  notation  is  precisely  Eq.  (4.24): 


G2(t,y;0,Y)  » “ | 


2^  ^ 2yY  's 


(D.6) 


where 


p(t)  » /t2+  (y-Y)  2 


(D.7a) 


0(t)  « /f27(y+y)2  . 


(D.7b) 


Continuing  to  follow  the  analysis  of  section  4 
we  find  our  two  forms  for  G: 


G(x,y;X,Y)  - G^  (x-X,y;0,Y)  + 2Gj^  (x,x+X;X,  0)  G2  (x+X,y;0 , Y) 
rx+X 

-2  (x,t;X,0)G,(t,y;0,Y)dt  (D.8) 

^x-X  t 
/■X+X 

-2  G (x,t;X,0)G,  (t,y;0,Y)dt.  (D.9) 

J x-X  1 


These  forms  are  identical  to  Eqs.  (4.7)  and  (4.8)  of 
the  memuscript.  In  Eq.  (D.8)  the  required  logarithmic 
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singularity  of  G is  apparent* in  the  first  term 

G (x-X,y;0,y).  As  a matter  of  fact,  it  was  proven  in 
2 

the  manuscript  (section  4b)  in  an  even  more  general 
form.  Combining  (D.5)  with  (D.8)  we  find  that 


G(x,0;X,Y)  H 0 


as  desired.  Also,  since  G^  and  G2  are  even  fvinc- 
tions  of  t,  the  integramd  of  (D.9)  is  odd  and  from 
(D.9)  we  see  that 


G(0,y;X,Y)  i 0. 


1 


Finally,  direct  differentiation  of  (D.8)  emd  (D.9) 
verifies  that  G is  indeed  a solution  of  the  homo- 
geneous  form  of  Eq.  (D.l)  when  (x,y)  (X,Y). 

Substitution  of  (D.3)  euid  (D.6)  into  (D.8)  yields 

G(x,y;X,Y)  • - W"  p 


2*  J » (x+x)  “ 


,x+X 

-2  — ^ P' 

2xX  n 


,1  /X^+X^-t^  . 
m^  2xX  ‘ 


. 1 f p .dt  (D.ii 

2^  Jp(t)  “ 2yY  s 

where  p emd  p are  as  defined  in  (D.7a)  and  (D.7b). 

The  situation  simplifies  somewhat  in  the  case  m«l»n. 
Then 


w 


54 


p (x-X) 


^ ^ ^ 1_  y^+Y^+ ix+K) 2 
2ir  2it  2yY 


0 (x^-X) 
0 (x+X) 


^x-X  27T  2yY  p(t) 


The  integral  may  be  written  in  a closed  form  consist- 
ing of  expressions  containing  powers  of  (x±X)  and 

2 2 

Cy±Y)  and  logarithms  of  (x±X)  + (y±Y)  [see 

reference  6,  pg.  205,  2.732]. 


I 


! 


\ 


55 


Appendix  S:  AN  ADDITION  FORMULA  FOR  RIEMANN  FUNCTIONS 

FOR  HYPERBOLIC  EQUATIONS 

As  mentioned  in  section  1 of  the  manuscript,  the 
motivation  for  this  thesis  in  elliptic  partial  differ- 
ential equations  came  from  recent  advances  in  the  theory 

of  Riemann  functions  for  hyperbolic  partial  differen- 
tial equations.  Because  of  the  nice  properties  of 
hyperbolic  equations  (e.g. , the  Riemann  function  is 
the  solution  to  a homogeneous  equation  and, thus, has 

no  singularity,  the  domain  of  influence  of  hyperbolic 
equations  is  bounded  by  two  real  cheuracteristics  and 

thus  finite  in  the  spatial  direction,  the  Riemamn 
function  does  not  depend  on  boundaur^  or  initial  con- 
ditions) it  is  much  sinpler  to  state  and  prove  an 
addition  formula  in  the  hyperbolic  case.  The  follow- 
ing theorem  and  corollary  may  be  found  in  reference 
3.  The  proof  consists  basically  of  a change  of  depen- 
dent variable  to  transform  the  equation  from  non-self- 
adjoint  form  into  self-adjoint  form  (as  was  done  in 
Appendix  C of  this  thesis)  and  direct  substitution 
and  differentiation  in  the  equation  (as  was  done  in 
section  3c  of  the  manuscript) . 

THEOREM.  If  Rj^(x,y,X,Y)  and  are 

the  respective  Riemeuin  functions  for 
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then  the  Riemann  function  R(x,y,X,Y)  for 

^xx“V  Ic3^(x)-C2(y)]U 

is  given  either  by 

R(x,y,X,y)  =»  Rj^(x,y-Y,X,0)  exp[|^  b2(t)dt] 
fx— X 

+ R, (x,t,X,0)  R,  (t,y,0,Y)  dt 

^ y-Y  ^ ^t 

or  by 

R(x,y,X,Y)  * R (x-X,y,0,Y)  exp [ [ b, (t)  dt] 

2 Jx 

ry-Y 

+ R2(t,y, 0,Y)  Rj^.  (x,t,X,0)  dt. 

^x-X  t 

COROLLARY.  Consider  the  eauation 


U + 4b, (r+s) [U  +U„]  + 4b. (r-s) [U  -U^] 
rs  i r s 2 r 3 


+ 4[Cj^(r+s)  - C2(r-s)]U  * 0. 


If  the  Riemann  function  of  (E.l)  is  Vj^(r,s,R,S)  when 

b2  ■ ■ 0 and  is  V2(r,s,R,S)  when 

then  the  Riemann  function  of  (E.l)  is  given  by 


V(r,s,R,S)  - V, (r,s,R,S)  exp  [|  b (t)  dt] 
^ J R+S  1 

rr-R+s-S  ^ ^ ^ 

j.  I tt  /r+s+t  r+s-t  R+S  R+Sx  . 

J r-R-s+S  1 ‘~2 2 r’ 


„ /t+r-s  t-r+s  R-S  -R+S  . , 

V2^(— 2— . — 2~  ' ~ '“T"  ) . 
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